Genetic sex validation for sample tracking in next-generation sequencing clinical testing.

OBJECTIVE: Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups. RESULTS: Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors (49.09%), samples from transgender participants (3.64%) and stem cell or bone marrow transplant patients (7.27%) along with undetermined sample mix-ups (40%) for which sample swaps occurred prior to arrival at genome centers, however the exact cause of the events at the sampling sites resulting in the mix-ups were not able to be determined.

Genetic Sex Validation for Sample Tracking in Clinical Testing.

Objective: Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups. Results: Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors, samples from transgender participants and stem cell or bone marrow transplant patients along with undetermined sample mix-ups.

Long read sequencing and expression studies of AHDC1 deletions in Xia-Gibbs syndrome reveal a novel genetic regulatory mechanism.

Xia-Gibbs syndrome (XGS; MIM# 615829) is a rare mendelian disorder characterized by Development Delay (DD), intellectual disability (ID), and hypotonia. Individuals with XGS typically harbor de novo protein-truncating mutations in the AT-Hook DNA binding motif containing 1 (AHDC1) gene, although some missense mutations can also cause XGS. Large de novo heterozygous deletions that encompass the AHDC1 gene have also been ascribed as diagnostic for the disorder, without substantial evidence to support their pathogenicity. We analyzed 19 individuals with large contiguous deletions involving AHDC1, along with other genes. One individual bore the smallest known contiguous AHDC1 deletion (∼350 Kb), encompassing eight other genes within chr1p36.11 (Feline Gardner-Rasheed, IFI6, FAM76A, STX12, PPP1R8, THEMIS2, RPA2, SMPDL3B) and terminating within the first intron of AHDC1. The breakpoint junctions and phase of the deletion were identified using both short and long read sequencing (Oxford Nanopore). Quantification of RNA expression patterns in whole blood revealed that AHDC1 exhibited a mono-allelic expression pattern with no deficiency in overall AHDC1 expression levels, in contrast to the other deleted genes, which exhibited a 50% reduction in mRNA expression. These results suggest that AHDC1 expression in this individual is compensated by a novel regulatory mechanism and advances understanding of mutational and regulatory mechanisms in neurodevelopmental disorders.

Proteogenomic Markers of Chemotherapy Resistance and Response in Triple-Negative Breast Cancer.

Microscaled proteogenomics was deployed to probe the molecular basis for differential response to neoadjuvant carboplatin and docetaxel combination chemotherapy for triple-negative breast cancer (TNBC). Proteomic analyses of pretreatment patient biopsies uniquely revealed metabolic pathways, including oxidative phosphorylation, adipogenesis, and fatty acid metabolism, that were associated with resistance. Both proteomics and transcriptomics revealed that sensitivity was marked by elevation of DNA repair, E2F targets, G2-M checkpoint, interferon-gamma signaling, and immune-checkpoint components. Proteogenomic analyses of somatic copy-number aberrations identified a resistance-associated 19q13.31-33 deletion where LIG1, POLD1, and XRCC1 are located. In orthogonal datasets, LIG1 (DNA ligase I) gene deletion and/or low mRNA expression levels were associated with lack of pathologic complete response, higher chromosomal instability index (CIN), and poor prognosis in TNBC, as well as carboplatin-selective resistance in TNBC preclinical models. Hemizygous loss of LIG1 was also associated with higher CIN and poor prognosis in other cancer types, demonstrating broader clinical implications. SIGNIFICANCE: Proteogenomic analysis of triple-negative breast tumors revealed a complex landscape of chemotherapy response associations, including a 19q13.31-33 somatic deletion encoding genes serving lagging-strand DNA synthesis (LIG1, POLD1, and XRCC1), that correlate with lack of pathologic response, carboplatin-selective resistance, and, in pan-cancer studies, poor prognosis and CIN. This article is highlighted in the In This Issue feature, p. 2483.

Clinical and molecular features of pediatric cancer patients with Lynch syndrome.

BACKGROUND: The association of childhood cancer with Lynch syndrome is not established compared with the significant pediatric cancer risk in recessive constitutional mismatch repair deficiency syndrome (CMMRD). PROCEDURE: We describe the clinical features, germline analysis, and tumor genomic profiling of patients with Lynch syndrome among patients enrolled in pediatric cancer genomic studies. RESULTS: There were six of 773 (0.8%) pediatric patients with solid tumors identified with Lynch syndrome, defined as a germline heterozygous pathogenic variant in one of the mismatch repair (MMR) genes (three with MSH6, two with MLH1, and one with MSH2). Tumor analysis demonstrated evidence for somatic second hits and/or increased tumor mutation burden in three of four patients with available tumor with potential implications for therapy and identification of at-risk family members. Only one patient met current guidelines for pediatric cancer genetics evaluation at the time of tumor diagnosis. CONCLUSION: Approximately 1% of children with cancer have Lynch syndrome, which is missed with current referral guidelines, suggesting the importance of adding MMR genes to tumor and hereditary pediatric cancer panels. Tumor analysis may provide the first suggestion of an underlying cancer predisposition syndrome and is useful in distinguishing between Lynch syndrome and CMMRD.

Implementation of preemptive DNA sequence-based pharmacogenomics testing across a large academic medical center: The Mayo-Baylor RIGHT 10K Study.

PURPOSE: The Mayo-Baylor RIGHT 10K Study enabled preemptive, sequence-based pharmacogenomics (PGx)-driven drug prescribing practices in routine clinical care within a large cohort. We also generated the tools and resources necessary for clinical PGx implementation and identified challenges that need to be overcome. Furthermore, we measured the frequency of both common genetic variation for which clinical guidelines already exist and rare variation that could be detected by DNA sequencing, rather than genotyping. METHODS: Targeted oligonucleotide-capture sequencing of 77 pharmacogenes was performed using DNA from 10,077 consented Mayo Clinic Biobank volunteers. The resulting predicted drug response-related phenotypes for 13 genes, including CYP2D6 and HLA, affecting 21 drug-gene pairs, were deposited preemptively in the Mayo electronic health record. RESULTS: For the 13 pharmacogenes of interest, the genomes of 79% of participants carried clinically actionable variants in 3 or more genes, and DNA sequencing identified an average of 3.3 additional conservatively predicted deleterious variants that would not have been evident using genotyping. CONCLUSION: Implementation of preemptive rather than reactive and sequence-based rather than genotype-based PGx prescribing revealed nearly universal patient applicability and required integrated institution-wide resources to fully realize individualized drug therapy and to show more efficient use of health care resources.

Transcriptional Reprogramming Differentiates Active from Inactive ESR1 Fusions in Endocrine Therapy-Refractory Metastatic Breast Cancer.

Genomic analysis has recently identified multiple ESR1 gene translocations in estrogen receptor alpha-positive (ERα+) metastatic breast cancer (MBC) that encode chimeric proteins whereby the ESR1 ligand binding domain (LBD) is replaced by C-terminal sequences from many different gene partners. Here we functionally screened 15 ESR1 fusions and identified 10 that promoted estradiol-independent cell growth, motility, invasion, epithelial-to-mesenchymal transition, and resistance to fulvestrant. RNA sequencing identified a gene expression pattern specific to functionally active ESR1 gene fusions that was subsequently reduced to a diagnostic 24-gene signature. This signature was further examined in 20 ERα+ patient-derived xenografts and in 55 ERα+ MBC samples. The 24-gene signature successfully identified cases harboring ESR1 gene fusions and also accurately diagnosed the presence of activating ESR1 LBD point mutations. Therefore, the 24-gene signature represents an efficient approach to screening samples for the presence of diverse somatic ESR1 mutations and translocations that drive endocrine treatment failure in MBC. SIGNIFICANCE: This study identifies a gene signature diagnostic for functional ESR1 fusions that drive poor outcome in advanced breast cancer, which could also help guide precision medicine approaches in patients harboring ESR1 mutations.

Genetic testing in ambulatory cardiology clinics reveals high rate of findings with clinical management implications.

PURPOSE: Cardiovascular disease (CVD) is the leading cause of death in adults in the United States, yet the benefits of genetic testing are not universally accepted. METHODS: We developed the "HeartCare" panel of genes associated with CVD, evaluating high-penetrance Mendelian conditions, coronary artery disease (CAD) polygenic risk, LPA gene polymorphisms, and specific pharmacogenetic (PGx) variants. We enrolled 709 individuals from cardiology clinics at Baylor College of Medicine, and samples were analyzed in a CAP/CLIA-certified laboratory. Results were returned to the ordering physician and uploaded to the electronic medical record. RESULTS: Notably, 32% of patients had a genetic finding with clinical management implications, even after excluding PGx results, including 9% who were molecularly diagnosed with a Mendelian condition. Among surveyed physicians, 84% reported medical management changes based on these results, including specialist referrals, cardiac tests, and medication changes. LPA polymorphisms and high polygenic risk of CAD were found in 20% and 9% of patients, respectively, leading to diet, lifestyle, and other changes. Warfarin and simvastatin pharmacogenetic variants were present in roughly half of the cohort. CONCLUSION: Our results support the use of genetic information in routine cardiovascular health management and provide a roadmap for accompanying research.

Sequencing of a central nervous system tumor demonstrates cancer transmission in an organ transplant.

Four organ transplant recipients from an organ donor diagnosed with anaplastic pleomorphic xanthoastrocytoma developed fatal malignancies for which the origin could not be confirmed by standard methods. We identified the somatic mutational profiles of the neoplasms using next-generation sequencing technologies and tracked the relationship between the different samples. The data were consistent with the presence of an aggressive clonal entity in the donor and the subsequent proliferation of descendent tumors in each recipient. Deleterious mutations in BRAF, PIK3CA, SDHC, DDR2, and FANCD2, and a chromosomal deletion spanning the CDKN2A/B genes, were shared between the recipients' lesions. In addition to demonstrating that DNA sequencing tracked a donor/recipient cancer transmission, this study established that the genetic profile of a donor tumor and its potential aggressive phenotype could have been determined before transplantation was considered. As the genetic correlates of tumor invasion and metastases become better known, adding genetic profiling by DNA sequencing to the data considered for transplant safety should be considered.

The Exceptional Responders Initiative: Feasibility of a National Cancer Institute Pilot Study.

BACKGROUND: Tumor molecular profiling from patients experiencing exceptional responses to systemic therapy may provide insights into cancer biology and improve treatment tailoring. This pilot study evaluates the feasibility of identifying exceptional responders retrospectively, obtaining pre-exceptional response treatment tumor tissues, and analyzing them with state-of-the-art molecular analysis tools to identify potential molecular explanations for responses. METHODS: Exceptional response was defined as partial (PR) or complete (CR) response to a systemic treatment with population PR or CR rate less than 10% or an unusually long response (eg, duration >3 times published median). Cases proposed by patients' clinicians were reviewed by clinical and translational experts. Tumor and normal tissue (if possible) were profiled with whole exome sequencing and, if possible, targeted deep sequencing, RNA sequencing, methylation arrays, and immunohistochemistry. Potential germline mutations were tracked for relevance to disease. RESULTS: Cases reflected a variety of tumors and standard and investigational treatments. Of 520 cases, 476 (91.5%) were accepted for further review, and 222 of 476 (46.6%) proposed cases met requirements as exceptional responders. Clinical data were obtained from 168 of 222 cases (75.7%). Tumor was provided from 130 of 168 cases (77.4%). Of 117 of the 130 (90.0%) cases with sufficient nucleic acids, 109 (93.2%) were successfully analyzed; 6 patients had potentially actionable germline mutations. CONCLUSION: Exceptional responses occur with standard and investigational treatment. Retrospective identification of exceptional responders, accessioning, and sequencing of pretreatment archived tissue is feasible. Data from molecular analyses of tumors, particularly when combining results from patients who received similar treatments, may elucidate molecular bases for exceptional responses.

Molecular Features of Cancers Exhibiting Exceptional Responses to Treatment.

A small fraction of cancer patients with advanced disease survive significantly longer than patients with clinically comparable tumors. Molecular mechanisms for exceptional responses to therapy have been identified by genomic analysis of tumor biopsies from individual patients. Here, we analyzed tumor biopsies from an unbiased cohort of 111 exceptional responder patients using multiple platforms to profile genetic and epigenetic aberrations as well as the tumor microenvironment. Integrative analysis uncovered plausible mechanisms for the therapeutic response in nearly a quarter of the patients. The mechanisms were assigned to four broad categories-DNA damage response, intracellular signaling, immune engagement, and genetic alterations characteristic of favorable prognosis-with many tumors falling into multiple categories. These analyses revealed synthetic lethal relationships that may be exploited therapeutically and rare genetic lesions that favor therapeutic success, while also providing a wealth of testable hypotheses regarding oncogenic mechanisms that may influence the response to cancer therapy.

Novel CYP2A6 diplotypes identified through next-generation sequencing are associated with in-vitro and in-vivo nicotine metabolism.

OBJECTIVES: Smoking patterns and cessation rates vary widely across smokers and can be influenced by variation in rates of nicotine metabolism [i.e. cytochrome P450 2A6 (CYP2A6), enzyme activity]. There is high heritability of CYP2A6-mediated nicotine metabolism (60-80%) owing to known and unidentified genetic variation in the CYP2A6 gene. We aimed to identify and characterize additional genetic variants at the CYP2A6 gene locus. METHODS: A new CYP2A6-specific sequencing method was used to investigate genetic variation in CYP2A6. Novel variants were characterized in a White human liver bank that has been extensively phenotyped for CYP2A6. Linkage and haplotype structure for the novel single nucleotide polymorphisms (SNPs) were assessed. The association between novel five-SNP diplotypes and nicotine metabolism rate was investigated. RESULTS: Seven high-frequency (minor allele frequencies ≥6%) noncoding SNPs were identified as important contributors to CYP2A6 phenotypes in a White human liver bank (rs57837628, rs7260629, rs7259706, rs150298687 (also denoted rs4803381), rs56113850, rs28399453, and rs8192733), accounting for two times more variation in in-vitro CYP2A6 activity relative to the four established functional CYP2A6 variants that are frequently tested in Whites (CYP2A6*2, *4, *9, and *12). Two pairs of novel SNPs were in high linkage disequilibrium, allowing us to establish five-SNP diplotypes that were associated with CYP2A6 enzyme activity (rate of nicotine metabolism) in-vitro in the liver bank and in-vivo among smokers. CONCLUSION: The novel five-SNP diplotype may be useful to incorporate into CYP2A6 genotype models for personalized prediction of nicotine metabolism rate, cessation success, and response to pharmacotherapies.

Postmortem genetic screening for the identification, verification, and reporting of genetic variants contributing to the sudden death of the young.

Each year in the United States, thousands of cases of sudden and unexpected deaths of infants, children, and young adults are assigned an undetermined cause of death after postmortem investigation and autopsy. Heritable genetic variants have been suggested as the cause of up to a third of sudden death (SD) cases. Elucidation of the genetic variants involved in SD cases is important to not only help establish cause and manner of death of these individuals, but to also aid in determining whether familial genetic testing should be considered. Previously, these types of postmortem screenings have not been a feasible option for most county medical examiners' and coroners' offices. We sequenced full exons of 64 genes associated with SD in the largest known cohort (351) of infant and young SD decedents using massively parallel sequencing at <$600 per sample. Genetic variants were assessed through literature review and clinical evaluation by a multidisciplinary consortium of experts. Thirteen individuals (3.7%), eight infants (2.8% of those <1 yr of age) and five children/young adults (7.0% of those >1 yr of age), were found to have a reportable genetic variant contributing to SD. These percentages represent an estimate lower than those previously reported. Overall yields and results likely vary between studies due to differences in evaluation techniques and reporting. Additionally, we recommend ongoing assessment of data, including nonreported novel variants, as technology and literature continually advance. This study demonstrates a strategy to implement molecular autopsies in medicolegal investigations of young SD decedents.

Genomic signatures of cooperation and conflict in the social amoeba.

Cooperative systems are susceptible to invasion by selfish individuals that profit from receiving the social benefits but fail to contribute. These so-called "cheaters" can have a fitness advantage in the laboratory, but it is unclear whether cheating provides an important selective advantage in nature. We used a population genomic approach to examine the history of genes involved in cheating behaviors in the social amoeba Dictyostelium discoideum, testing whether these genes experience rapid evolutionary change as a result of conflict over spore-stalk fate. Candidate genes and surrounding regions showed elevated polymorphism, unusual patterns of linkage disequilibrium, and lower levels of population differentiation, but they did not show greater between-species divergence. The signatures were most consistent with frequency-dependent selection acting to maintain multiple alleles, suggesting that conflict may lead to stalemate rather than an escalating arms race. Our results reveal the evolutionary dynamics of cooperation and cheating and underscore how sequence-based approaches can be used to elucidate the history of conflicts that are difficult to observe directly.

A massive expansion of effector genes underlies gall-formation in the wheat pest Mayetiola destructor.

Gall-forming arthropods are highly specialized herbivores that, in combination with their hosts, produce extended phenotypes with unique morphologies [1]. Many are economically important, and others have improved our understanding of ecology and adaptive radiation [2]. However, the mechanisms that these arthropods use to induce plant galls are poorly understood. We sequenced the genome of the Hessian fly (Mayetiola destructor; Diptera: Cecidomyiidae), a plant parasitic gall midge and a pest of wheat (Triticum spp.), with the aim of identifying genic modifications that contribute to its plant-parasitic lifestyle. Among several adaptive modifications, we discovered an expansive reservoir of potential effector proteins. Nearly 5% of the 20,163 predicted gene models matched putative effector gene transcripts present in the M. destructor larval salivary gland. Another 466 putative effectors were discovered among the genes that have no sequence similarities in other organisms. The largest known arthropod gene family (family SSGP-71) was also discovered within the effector reservoir. SSGP-71 proteins lack sequence homologies to other proteins, but their structures resemble both ubiquitin E3 ligases in plants and E3-ligase-mimicking effectors in plant pathogenic bacteria. SSGP-71 proteins and wheat Skp proteins interact in vivo. Mutations in different SSGP-71 genes avoid the effector-triggered immunity that is directed by the wheat resistance genes H6 and H9. Results point to effectors as the agents responsible for arthropod-induced plant gall formation.